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Abstract 

We perform a perturbative calculation of the physical observables, in particular pseudo- 
Hermitian position and momentum operators, the equivalent Hermitian Hamiltonian op- 
erator, and the classical Hamiltonian for the 'PT-symmetric cubic anharmonic oscillator, 
H = + ^fi 2 x 2 + iex 3 . Ignoring terms of order e 4 and higher, we show that this system 
describes an ordinary quartic anharmonic oscillator with a position-dependent mass and 
real and positive coupling constants. This observation elucidates the classical origin of the 
reality and positivity of the energy spectrum. We also discuss the quantum-classical cor- 
respondence for this PT-symmetric system, compute the associated conserved probability 
density, and comment on the issue of factor-ordering in the pseudo-Hermitian canonical 
quantization of the underlying classical system. 

PACS numbers: 03.65.-w 
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1 Introduction 



PT-symmetric quantum mechanics was originated by the observation, initially made by Bessis 
and Zinn- Justin, that the Hamiltonian for a cubic anharmonic oscillator: 

with /i, e G R has a real, positive, and discrete spectrum. During the past six years there have 
appeared a number of publications PQ-[E] exploring the properties of the Hamiltonian (JJJ. Yet 
the nature of the physical system described by this Hamiltonian has not been clarified. The 
present article aims at addressing this basic issue. We will achieve this aim by computing the 
physical observables, the localization probability density, and the underlying classical Hamil- 
tonian for this system. This is the first example of a PT-symmetric quantum system with 
configuration space R that allows for such a computation. 

As the main technical tools used in our analysis have been developed in the study of the 
spectral properties of PT-symmetric Hamiltonians, we will include here a brief account of the 
relevant developments. 

The first convincing numerical evidence supporting the reality and positivity of the spectrum 
of (JTJ was provided by Bender and Boettcher £Q who made the important observation that this 
Hamiltonian was PT-symmetric. Among dozens of publications on the subject that followed 

was the article by Dorey, Dunning, and Tateo jU] that provided the first mathematically 
rigorous proof of the spectral properties conjectured by Bessis and Zinn- Justin (See also [TU]). 
From a physicist's point of view, a more important development was the idea, put forward by 
Bender and his collaborators PQE]; that such PT-symmetric Hamiltonians might be used as 
the Hamiltonian operator for an extended/generalized quantum theory. 

The main obstacle for realizing this idea was that a non-Hermitian Hamiltonian such as 
(PQ) generated a nonunitary time-evolution. This was not compatible with the conventional 
probabilistic interpretation of quantum mechanics. The resolution of this problem came as a 
by-product of the attempts to characterize the non-Hermitian operators having a real spectrum, 

mm 

Ref. [13J lists the necessary and sufficient conditions that ensure the reality of the spectrum 
of a diagonalizable operator. 1 Among these is the condition that H must be Hermitian with 
respect to a positive-definite inner product (•,•)+. This inner product is generally different from 
the defining inner product (-, •) of the (reference) Hilbert space 7i in which the operator H acts. 
It can be conveniently expressed in terms of a positive-definite (metric) operator r/ + : 7i — > 7i 
according to [TJ] 

(;•)+= (;V+.). (2) 

The condition that H be Hermitian with respect to (-, ■)+, i.e., (•, H-) + = (H-, •)+, is equivalent 
to 7/_|_-pseudo-Hermiticity of H, [T2]. This means that rj + belongs to the set Uh of all Hermitian 

view of the requirements of the standard quantum measurement theory physical observables in general 
and the Hamiltonian in particular must be diagonalizable operators |17| . 
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invertible operators n : 7i — > 7i satisfying 2 ^2] 

H ] = r]Hr]- 1 . (3) 

An interesting property of the set Uh of all metric operators r\ is that to each pair (771,772) of 
elements of Uh there corresponds a symmetry generator r]2 1 i]i of H, ^2]- Furthermore VT- 
symmetric Hamiltonians H that act in TL = L 2 (R), e.g., (fl|). are P-pseudo- Hermitian, i.e., 
P G This in turn implies that if H has a real spectrum then V~ 1 i] + = Vi] + commutes 
with H. The construction of the physical Hilbert space 7i p h y s that is based on the CPT-inner 
product makes an implicit use of this observation. As shown in ^Bj for theories defined 
on R and more recently generalized in ^H] to theories defined on a complex contour, the C 
operator introduced in [IR is related to the metric operator r/ + according to 

C = Vr 1+ = rj^V, (4) 

and the CPT-inner product is precisely (•,•)+■ 

The recent approximate calculations of C for the anharmonic oscillator and its analogs [HI 
20 have also revealed the practical significance of the factorization of C. These calculations 
are based on Eq. (j3J and the observation that (being a positive-definite operator) rj + admits 
an exponential representation, 

V+ = e- Q (5) 

where Q is a Hermitian operator. 

The metric operator r)+ (which is generally unique up to symmetries of H, [13 HEl EI]), 
not only determines the structure of the physical Hilbert space but it fixes the observables of 
the theory as well, [22 123 EI- By definition (Def. 1.) Physical observables are the Hermitian 
operators acting in the physical Hilbert space TC p hy S , {22H2HHJZ!, i-e., A : H — ► TC is an observables 
if 3 

(;A.) + =(A;.) + . (6) 

Alternatively, Physical observables A are n + -pseudo-Hermitian operators acting in 7i, i.e., A' = 
V+AV+ 1 - 

In order to see the central role played by the metric operator rj + in the construction of the 
observables, we recall that as an operator mapping H + onto TC the unique positive square root 
p = ■ s /rj, of r] + is a unitary operator [2EllIIj, i-e., 

(P ■,/>•) = (•, •>+■ (7) 

2 Here and throughout this paper the adjoint of an operator A : Hi — * H.2 between two Hilbert spaces Hi 
(with inner product (•, and H2 (with inner product (•, ^2) is defined to be the unique operator : H2 — > Hi 
satisfying (;A-) 2 = (^ f -,-)i. 

3 As argued in [221 identifying observables with CPT-invariant operators as initially done in 1|| leads to a 
dynamical inconsistency. The latter is avoided if one modifies this definition as proposed in [33] . This modified 
definition is equivalent to Def. 1 above for symmetric Hamiltonians H (satisfying (x\H\x'} = (x'\H\x)) and 
cannot be applied for nonsymmetric Hamiltonians |25| . 
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Hence the Hermitian operators O acting in ri p ^ ys (i-e., the physical observables) may be obtained 
from the Hermitian operators o acting in TC according to 

= p- 1 op. (8) 

This is also consistent with the condition that H is related to a Hermitian operator h : 
TC — > TC by a similarity transformation, 

h = pHp-\ (9) 

The mapping p : TC p h ys — ► W establishes the unitary equivalence of the PT-symmetric quantum 
system S PT having 7i p h ys , H, and O as the physical Hilbert space, the Hamiltonian, and the 
physical observables and the quantum system S having TC, h, and o as the physical Hilbert 
space, the Hamiltonian, and the physical observables, respectively, |22j [T^J . Spt and S describe 
the same physical system because the physical quantities such as the expectation values and 
transition amplitudes are independent of the choice of Spt and S. 

The advantage of the PT-symmetric description provided by Spt over the Hermitian de- 
scription provided by S is that unlike H, the Hermitian Hamiltonian h is generally nonlocal. 
This advantage is however balanced by the disadvantage that the physical (pseudo-Hermitian) 
position X and momentum operators P of Spp are also generally nonlocal. These operators 
are defined by [TH] 

X:=p- X xp, P:=p~ l pp, (10) 

where x and p are the conventional position and momentum operators. The main advantage 
of the Hermitian description is that it provides means for identifying the underlying classical 
system, [T7|. The classical Hamiltonian is obtained by expressing h in terms of x and p and 
replacing the latter with the classical (real-valued) position x c and momentum p c observables. 
In general this yields an expression that may involve powers of h. The classical Hamiltonian 
H c is then obtained by evaluating this expression in the limit H — *■ 0, i.e., assuming that this 
limit exists, 

H c (x c ,p c ) := \imh(x c ,p c ). (11) 

The initial Hamiltonian H may be recovered by performing the so-called ?7 + -pseudo-Hermitian 
canonical quantization of H c and adopting an appropriate factor-ordering prescription [TT] . 
Disregarding the complications due to the factor-ordering problem and assuming that H c is an 
analytic function of x c and p c , we have 

H C (X,P) = h(X,P) = h(p~ l x p, p- 1 p p) = p- 1 h(x,p)p = H. (12) 

Having introduced the ?7 + -pseudo-Hermitian position operator X we can also address the 
issue of determining the conserved probability density p for the localization of the system in 
the configuration space. This requires the identification of the physical localized states of the 
system. Being (the generalized j2Zj) eigenvectors of X, the localized state vectors are given by 

\& ) )=p- 1 \x), (13) 
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where \x) are the conventional position eigenvectors. The conserved probability density asso- 
ciated with a given state vector ip e 7i p h ys has the form [T7] 

e [x) = N^l^ix)] 2 , (14) 

where ty(x) is the physical position wave function for the state vector ip, i.e., 

*(x):=$ x \,J>) + =(x\pil>), (15) 

(•|-) is the usual L 2 -inner product on Ti = L 2 (R), and N := {tp,ip) + = \^{x)\ 2 dx. A 

An important feature of the exponential representation (J5J) of the metric operator 77+ is that 
it reduces the calculation of p and p _1 to that of Q, for 

p ±l = e^ 2 . (16) 

We will use this observation together with the approach pursued in jH] to perform a perturbative 
calculation of X, P, h, H c , and g for the PT-symmetric Hamiltonian ([T|). 

Because we are interested in the issue of finding the classical limit of the PT-symmetric 
theory based on the Hamiltonian (0) we wish to retain the factors of h. However, for the sim- 
plicity of the calculations and ease of the comparison with the known results, we will introduce 
and employ the following dimensionless quantities. 



x := £ 1 x 



p := ih-'p, (17) 
ftr l \pmu, e:=£ 5 h~ 2 me, (18) 



M := I n ym p 
H := ^p 2 + ^x 2 , H i: =*x 3 , (19) 

1 M 2 

H := £ 2 h- 2 mH = -p 2 + —x 2 + iex 3 = B. + eB. u (20) 
where I is an arbitrary length scale which may be taken as p 2 /e. Clearly, we have [x, p] = i. 



2 Calculation of Q 

In jH] the authors outline a perturbative calculation of Q for the Hamiltonian (|2U|) taking e as 
the perturbation parameter. They use the identities [C, VT] = and C 2 = 1 to infer that as a 
function of x and p, Q must be even in x and odd in p. Furthermore, imposing [C, H] = and 
making use of the fact that Hi is an imaginary cubic potential, they find the operator equation 

2ee Q E 1 = [e Q 1 E} 1 (21) 

4 The physical position wave functions evolve in time according to the Schrodinger equation with h being the 
Hamiltonian operator. In general, one can express h in the form p 2 /(2m) + W where W is a nonlocal potential 
(an infinite series in p with a;-dependent coefficients). This in turn implies that the probability current density 
that together with g satisfies the continuity equation has a nonlocal dependence on ^(x); it is not given by the 
standard formula, unless H is Hermitian. This is especially significant in the study of tunnelling and scattering 
for pseudo-Hermitian Hamiltonians (having scattering states). 
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and that Q may be expanded in an odd power series in e, 

Q = Q 1 e + Q 3 e 3 + Q 5 e 5 + ■■■ , (22) 

where Q2i+i = Q2i+i(x,p) with i = 0, 1,2, ••• are e-independent. Next, they expand e Q in 
power series in e, substitute the result in (jUJ), and demand that this equation be satisfied at 
each order of e. This yields a series of operator equations that they iteratively solve for C^i+i- 
The operator equations whose solution yield Q21+1 may be more conveniently obtained from 
the ?7 + -pseudo-Hermiticity of H, 

Ht^+Hr^ 1 . (23) 
Substituting r/ + = e~® in this equation and noting that W = H — eHi, we have 

H -e" Q H e Q = e(H 1 + e- Q H 1 e Q ). (24) 

Next, we employ the Baker-Campbell-Hausdorff identity, 

e- A B e A = B+ [B, A] + A], A] + A], A], A] + ■ ■ ■ , (25) 

(where A and B are linear operators), to express (j2*4^) as 

- [Ho, Q] ~ ^ [[H , Q],Q}-^ [[[H , Q],Q],Q] = e f 2H X + [Hi, Q] + i [\E lf Q],Q]+ 

I[[[ Hl ,Q],Q],Q] + ---V (26) 

Now, in view of (|22p. we can easily identify the terms in (f2*U|) that are of the same order in 
powers of e. Enforcing (|2*S|l at each order, we find the desired operator equations for Qzi+i- 
Matching the terms of order e, e 2 , ■ ■ ■ , e 5 , we find in this way the following independent operator 
equations which agree with those obtained in jS]. 5 

[H ,Qi] = -2H 1; (27) 

[H ,Q 3 ] = -lm 1 ,Q 1 ],Q 1 ], (28) 
o 

[H ,Q 5 ] = ~\ ([[Hi,Qi],Q3] + [[H 1 ,Q 3 ],Qi]) + -i-[[[[H 1 ,Q 1 ],Q 1 ],Q 1 ],Q 1 ]. (29) 

o •)()() 

The higher order terms in e similarly yield operator equations for Q21+1 with % > 3. As noted 
in jH], one can iteratively solve these equations to obtain Q2i+i- 
A variation of the approach of jH] is to substitute the ansatz 6 

i+l 

Q 2 m= ]Tc^{x 2 ^ 2fc+1 } (30) 

j,fc=0 



5 These equation are obtained at the orders e,e 3 and e 5 , respectively. 
6 Here {•, •} stands for the anticommutator, {A, B} = AB + BA. 
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in the operator equations for Q2i+i and to solve for the coefficients c^. In this way we have 
found the following solutions for (j2Z|) and (J2BJ), respectively. 



Qi 
Q 3 



l 

4 



M 



10 



4 
3 

32 



p 3 + M 2 {x 2 ,p} 
5 



4 



p + 2m x p x 



(31) 



, r P 5 + o ^ 2 { x2 ; P 3 } + >i 4 {x 4 , p} + 2.M 2 p 
15 3 



128 



P + 



40 



xp x 



2 2 
X p X 



32 



(32) 



15.M 10 3m 8 " M v " .M c 

These confirm the results of jS] except for the coefficient of the last term in (J32|) . We have 
checked the validity of (J32|) by inserting this relation in (J28|) and affecting both sides of the 
resulting equation on the function /i(x) = x. Using the fact that in the x-representation 
p = we could easily perform the necessary calculations (without having to use any 

commutation relations) and checked the validity (|32j) . 

In fact, we can obtain the coefficients Cijk using this method. In order to do this we can 
substitute (|3T)j) in the operator equations for Q2i+i, (rather than trying to use the complicated 
commutation relations for powers of x and p) affect both sides of these equations on / n (x) = x n , 
and demand that they are equal for all n — 0, 1, 2, 3, • • • . 



3 The Equivalent Hermitian Hamiltonian 

Having obtained Q, we can easily calculate the Hermitian Hamiltonian 

h = pEp- 1 

associated with the dimensionless Hamiltonian H. Using (jlfij) . (|25jl . and ()33)1 . we have 

[E,Q],Q],Q}+--- . (34) 



f33) 



h = H + -[H,Q]-, 2 , 22 



' [[H,Q],Q]+ ' 



3!2 3 ' 



Now, in view of (j2*U|) and (|22p. it is very easy to identify the perturbative expansion of h, i.e. 
find e-independent operators such that 



This yields 
h (o 

h< 2 

h (3 

h( 4 
h (5 



i=0 



h« = H 1 + -[H ,Q 1 ], 
~[Hi,Qi] + lmo,Qi},Qi], 

- [H , Q 3 ] + ~ [[H 1; Q 1 ], Qi] + 1 [[[H , Qi], Qi], Qi], 

\\ni,Q*)-^[[\RiMMQil 

\ [Ho, g 5 ] + 1 ([[Hi, Qx], q 3 ] + [[h 1; g 3 ], Qi]) + ^ 



H 
1 

2 
1 



(35) 

(36) 
(37) 
(38) 
(39) 

[H ,Q 3 ],Qi],Qi]- (40) 
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In view of the fact that Qi, Q3, Q5 are Hermitian while Hi is anti-Hermitian, it is not difficult 
to see that the terms contributing to hy' with even j are Hermitian while those contributing to 
hv') with odd j are anti-Hermitian. The fact that h is a Hermitian operator then suggests that 
the hy> with odd j must vanish. There is another argument supporting this expectation namely 
that because Hi is a cubic potential, the perturbation series for the ground state energy of H 
and consequently (the isospectral operator) h must only include even powers of the perturbation 
parameter e, [3]. 

Using (|27p. (|2Sj) . ([29)1 . we can easily show that indeed hW, h^, and h^ 5 ^ vanish identically. 
This may be viewed as a consistency check of our calculations. The perturbative expansion of 
h valid up to and including terms of order e 5 is, therefore, given by 

h = H + h(V + h(V + O(e 6 ), (41) 
h( 2 ) = l[Hi,Qi], (42) 

h« = \ [Hi,g 3 ] - ^ [[[Hi,Qi],Qi],Qi], (43) 

where we have made use of (|27j). Next, we use (|3*Tjl and (|3*2*Jl to obtain the explicit form of h^ 2 ^ 
and hy>. After a lengthy calculation, we find 

[Hi,Qi] = ^ ({x 2 ,P 2 } + ^ 4 + ; (44) 
[#1, Qa] = — \k (16{x 2 , p 4 } + 15.m 2 {x 4 , p 2 } + 64 p 2 + 6m 4 x 6 + 76m 2 x 2 ) , (45) 

M 

[[[Hi, Qi],Qi], Qi] = (8p 6 - 8m 2 {x 2 , p 4 } + 9m 4 {x\ p 2 } - 68a4 2 P 2 + 10^( 6 x 6 + 28m 4 x 2 ) . 

(46) 



Therefore, in view of (|42j) and (|43|) . 

h (2) = ^-J{*\p 2 } + M 2 x* + l), (47) 



2m a y 1 3 

h^ = A(V-W{x 2 ,p 4 }-^.M 4 {x 4 ,p 2 }-^ ^ 

A simple application of (|41|) is in the calculation of the energy eigenvalues E n of the Hamilto- 
nian H. If we denote by \n) the normalized eigenvectors of the harmonic oscillator Hamiltonian 
H , then we can easily calculate E„ up to and including terms of order e 3 . This is done using 
the first order Rayleigh-Schrodinger perturbation theory which yields 

E n = M (n + h + (n\h^\n) + 0(e 4 ). (49) 

Substituting (|47|) in this relation and doing the necessary algebra, we find 

E n = M {n + ~) + -i- (30n 2 + 30n + 11) e 2 + C(e 4 ). (50) 
2 oM 
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This is in complete agreement with the earlier calculations reported in |H[7j. 

Next, we use (J 17)) - (|19|) to obtain the expression for the unsealed Hermitian operator h that 
is associated with the original Hamiltonian H . This results in 

7 P 2 . 1 2 2 i ^ ( 1 r 2 2t 24, 2/ i l 2 \ 2 . 2 /j> 6 9,U 2 r 2 4i 



2m 2 2\r \m dm J fi lz \m d m z 

51/i 4 r 4 21 7/i 6 6 2 69^/i 4 2^ 4 ^/ 6N 

-{x 4 ,p 2 } ^-x 6 £-p 2 — x 2 e 4 + 0(e 6 ) (51) 



8m 4 2m 2 2m 

+ i^V + — f {* 2 ,p 2 } +px 2 p+ 6 2 + — (t- - -^{x 2 ,p 4 } 

2m 2 m/i 4 \ ' 2 J pP \m 3 16m 2 ' 

81/^ 2 2 2 33/i. r 4 2l ^9/A 2 2 2 ^A 4 &\ 4 , /o/ 6n /rn\ 

■p x p — — — {x ,p ) — x p x — x I e + C(e ), (52) 



8m 2 ' 16m ' 8m 
where we have used the identities 

2 ^ r 2 2i f-2 222 ^ r 4 2i a fc2 2 222 ^ r 2 4i 4 fc2 2 

pi p — -{x ,p } = a , xpx— —\x ,p } = 4ft x , p x p — -{x ,p } = An p . 

Note that if one does not truncate the perturbation expansion of h, one finds that it is an 
infinite series in powers of p. This confirms the assertion that the Hermitian Hamiltonian for a 
non-Hermitian Hamiltonian with a real spectrum is in general a nonlocal (pseudo-differential) 
operator, f2~E\ IT7]. A remarkable property of the cubic anharmonic oscillator (0), is that the 
corresponding Hermitian Hamiltonian h turns out to be a local (differential) operator once one 
truncates its perturbation expansion. This is not generally the case. 



4 Physical Observables 

The calculation of the physical observables O : 7i p h y s ^ P h ys mimics that of h. As we discussed 
in Section 1, because the reference Hilbert space 7i for the system is L 2 (IR), the observables O 
are obtained from the Hermitian operators o : L 2 (R) — > L 2 (M) according to JHJ). Substituting 
(dnj) in this relation and using (j2SJ), we have 

O = o - \[o, Q] + -L[[o, Q), Q] - Q),Q),Q)±.... (53) 

Moreover due to the particular e-dependence of Q as given by we can easily determine 
the following perturbation expansion for O. 

= o-±[o,Q 1 ]e+±[[o,Q 1 ],Q 1 ]e?-± f[o, Q 3 ] + ^[[[o, Q 1 ], Q 1 ], Qj) e 3 + 0(e'). (54) 

Next, we calculate the dimensionless r] + -pseudo- Hermitian position and momentum opera- 
tors, 

X:= p~ 1 xp = r 1 X, P := p- 1 pp = £h- 1 P. (55) 
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This is done by substituting x and p for o in (|54|). Doing the necessary calculations, we obtain 

X = x + ^ ( P 2 + \ m 2 A e + — 6 ({x, p 2 } - m 2 ^) e 2 + 0{e% (56) 

P = P --L{ x ,p} £ + -L (2p 3 -i M 2 {^, V }je 2 + 0(E 3 ). (57) 

We can directly read the expression for the ?7 + -pseudo-Hermitian position operator X and 
momentum operator P from these equations provided that we let X — * X, P — ► Pj sjm, x — » x, 
p — > p/ y'm, .m — > fx, e — > 6. As expected, X and P do not involve h. 

Equations (JHH)l and (|57j) show that, as operators acting in L 2 (R), X and P are not Hermitian. 
The fact that by construction they are ?y + -pseudo-Hermitian implies that as operators acting 
in 7iph ys they are Hermitian, [T2j. Furthermore, these operators furnish an irreducible unitary 
representation of the Heisenberg-Weyl algebra, [X, P] = ih, on the physical Hilbert space 
7iph ys - They form an irreducible set of basic operators for the quantum system, i.e., other 
observables may be constructed as power series in X and P. For instance, we can express the 
Hamiltonian JIJ according to 

H = — + -fi 2 X 2 + 1 ( [X 2 P 2 } + PX 2 P + ^ m ^ 2 x 4 ^] e 2 + — ( — - ^^—{X 2 P 4 } 
2m 2 m/i 4 \ ' 2 J /i 12 \m 3 16m 2 

p2x2p2 ~ I^ {x4 ' p2} ~ x2p2x2 ~ IT x ") e4 + ° (e6) ' (58) 

where we have made use of 0, (|10|). and (|52|1. This is the manifestly Hermitian representation 
of the original Hamiltonian (|T)). 

Another interesting implication of Eqs. (jSE)) and (JH7|) is that if e 7^ 0, the physical position 
(X) and momentum (P) operators do not satisfy the transformation rules of the usual position 
(2) and momentum (p) operators under V and T separately, 

VXV ^ -X, PPP ^ -P, TXT ^ T, TPT ^ -P. 

However, they share the same transformation rule under VT , 

VT XVT = -X, VTPVT = P. 

This is consistent with the fact that unlike V and T, VT is an antilinear r/ + -pseudo-unitary 
operator [2Ij- 7 In particular, it implies that, as an operator acting in 7i p hy S , VT is an antilinear 
unitary operator. This in turn implies, in view of Wigner's classification of symmetries in 
quantum mechanics [2E1, that unlike V and T, VT defines a physical symmetry of the quantum 
system. The fact that V does not correspond to a physical symmetry was to be expected, for 
its definition is intertwined with that of x which is not a physical observable for the system 
unless e = 0. 

7 This can be easily checked using the approach of • 
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5 The Classical Limit 



The phase space of the underlying classical Hamiltonian for the cubic anharmonic oscillator 
(JTJ is clearly M?. Having calculated the Hermitian operator h we can determine the classical 
Hamiltonian H c for this system using (fTTJ) . In view of (JHT|) , the evaluation of the limit in (fTT]) 
is trivial. Up to and including terms of order e 5 , H c is given by 

Pc i 1 ,,2^2 i ^ / 2 2 2 2 , J \ .2 



2 fpl 18 M 2 ,4_^42_M,j + C )f,- ( r,(|) 



^12 y m 3 m 2 c c c c ^ 

We shall first explore the consequences of neglecting the terms of order e 4 and higher. Then 
we can express H c in the form 

He = 777#^ + ^c + S^ + ^ 4 ), (60) 



777 

M(x c ) := 42 - =m(l-3/i-V^) + 0(6 4 ). (61) 

Therefore, for sufficiently small e, if c describes the dynamics of a point particle with a position- 
dependent mass M(x c ) that interacts with a quartic anharmonic potential. This statement 
provides a physical interpretation of the original PT-symmetric cubic anharmonic oscillator 
((TJ. Obviously, this is a valid approximation as long as we can neglect the contribution from 
the terms of order e 4 and higher, 0(e 4 ) ~ 0. 

Under this assumption, H c takes non-negative values; the classically allowed energies E 
are non-negative. This is the classical analog of the fact that the PT-symmetric quantum 
Hamiltonian ((TJ) has a positive spectrum. Moreover, it is not difficult to show that the classical 
orbits in the phase space for the Hamiltonian (|50j) are ellipses determined by 



2m ' V 2 ' yu 4 / 

The coupling of the energy E and the perturbation parameter e is an indication that the above 
approximation is valid for low energies, i.e., 

:=^V 2 - (63) 
o 

The inclusion of the terms of order e 4 distorts the above picture. However, as long as 
condition (jHSJ) holds the classical (phase-space) orbits are closed curves. Figure [T] shows the 
graph of such orbits. 

If we perform the ?7 + -pseudo- Hermitian quantization of the classical Hamiltonian ([59)1 . 
namely let x c — > X, p c — > P, and {•, -} c — > — z/i _1 [-,-], where {•, -} c is the classical Poisson 
bracket, we recover the expression (}58|) for the original Hamiltonian (fljl. provided that we 
adopt the correct factor-ordering prescription. This observation underlines the importance of 
the issue of factor-ordering ambiguity in pseudo-Hermitian and in particular "PT-symmetric 
quantum mechanics. 
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Figure 1: Graph of the orbits in phase space for the Hamiltonian (|59jl with 
C(e 6 ) neglected, m = \i = 1, e = 0.1, and E = 1 (dashed-dotted curve), 
E = 5 (dashed curve), E = 8 (dotted curve), and i? = 10 (full curve). The 
horizontal and vertical axes are those of x c and p c , respectively. Note that 
E* = 50/3 ~ 16.7. Hence our perturbative calculation of the classical orbit 
for E = 10 is not as reliable as that for E — 1 and E = 5. In particular the 
elliptic shape of the E = 1 orbit is consistent with Eq. ([62)1 . 

6 The Conserved Probability Density 

The expression ()14j) for invariant probability density for the localization of the quantum system 
under consideration involves the physical wave function Given a state vector if) e Ti. p hy S , 

the perturbation expansion for the corresponding physical wave function is given by 



me 



-9/2 1 



fc=0 



;-i) fc Q fc 

2 k k\ 



iP(x) - ^(xlQ^e + Ux\Qj\ 



\ «*| W> + k kAQWf)) ^ + O(e'), (64) 



where we have used (I16|) and ()22|) . We can obtain the explicit form of the terms appearing on 



the right-hand side of JSK using (|3TJ), (J32J), (fTTfl - (fTTHl. and the identity 
result may be expressed as 



ih^-(x\. The 



= (1 + eLi + e 2 L 2 + e 3 L 3 )^(x) + C(e 4 



(65) 



where 



-501, 



3 •= 



o - 4 * 



32ft 4 d 5 Wh 2 n 2 
1 L_ 

15m 2 dx 5 3m 



2 ^ 3 ~ 48 ^ :_ 71 



2H 2 d 3 + 2 2 d + 
3m dx 3 dx 



d 3 



2 " , o™ « \ o..4 ( „4 « , o™3 



x z —, + 3x—r) - 2^ [ x 4 — + 2x 

dx a aar / \ dx 



8ftV d_ 
m da; 
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Having obtained the general form of the physical wave function we can calculate the invariant 
probability density g according to (fT4|) . Figures 3 and 4 show the plots of g for ip(x) = e~ x ' 2//2 
and xe~ x2 l 2 . 



7 Conclusion 

We have performed a comprehensive study of the physical content of the PT-symmetric quan- 
tum system based on the non-Hermitian cubic anharmonic oscillator We showed how the 
general ideas developed within the framework of the pseudo-Hermitian quantum mechanics 
may be applied to this model. The result is an explicit characterization of the corresponding 
Hermitian Hamiltonian, physical observables, probability density, and the underlying classi- 
cal system. The only other PT-symmetric system (with an infinite-dimensional state space) 
for which a similar treatment has been possible is the PT-symmetric square well ^7]. An 
important difference between the latter system and the anharmonic oscillator ((TJ is that the 
effects of non-Hermiticity of this oscillator do survive the classical limit; non-Hermiticity is not 
a by-product of the (pseudo-Hermitian) quantization. 

Neglecting forth and higher order terms in our perturbative treatment, we showed that the 
PT-symmetric cubic anharmonic oscillator ([!} describes a point particle having a position- 
dependent mass and interacting with a real quartic anharmonic potential. This provides a 
classical justification for the positivity of the spectrum of (JTJ). The same argument applies to 
the cases where we should keep the terms of order up to (and including) five. 

The pseudo-Hermitian quantization of the classical Hamiltonian defined by the appropriate 
metric operator together with a particular factor-ordering prescription yields the original local 
PT-symmetric Hamiltonian while the usual canonical quantization of the same classical Hamil- 
tonian with the appropriate factor-ordering prescription leads to the corresponding equivalent 
nonlocal Hermitian Hamiltonian. 

The approach pursued in this paper may be applied to other PT-symmetric and non- 
PT-symmetric non-Hermitian Hamiltonians with a real spectrum. In general, however, the 
nonlocality of the corresponding equivalent Hermitian Hamiltonian may manifest itself at each 
order of the perturbation theory. This has already been the case for the PT-symmetric square 
well studied in ^7j. In view of the results of 8J, the same is the case for the PT-symmetric 
cubic potential, i.e., (JTJ) with /z = 0. An interesting subject of future study is to extend the 
approach pursued here to the field theoretical analog of (0). Such a study should reveal the 
structure of the underlying classical field theory. 



Note: After the online announcement of the preprint of this article ( quant-ph /0411137[) 



Hugh Jones sent me his preprint: quant-ph/0411171 in part of which he also studies the PT- 
symmetric cubic anharmonic oscillator. 
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Figure 2: Graph of the invariant probability density g for ip(x) = 
e~ x2 l 2 , H — m — n — 1 and e = (full curve) e = 0.2 (dashed 
curve), and e = 0.25 (dotted curve). 
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Figure 3: Graph of the invariant probability density g for ip(x) = 
xe -x 2 /2^ fa _ m _ ^ _ i anc j £ _ g curve ) e — o.l (dashed 

curve), and e = 0.15 (dotted curve). 
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Erratum 

There is a factor of 2 error in Eq. (61) which was noticed after the publication of the paper. 
Correcting this error leads to minor changes in Eqs. (62) and (63). The corrected equations are 

M (*c) := TTBTFW = m(l ~ 6/i" 4 e 2 x 2 ) + O(e'), (61) 
fi+(4+ C ^H-§^ = £, (62) 
^/i 6 e- 2 . (63) 

Eq. (62) shows that the distortion of the elliptic shape of the phase space orbits of the un- 
perturbed (harmonic oscillator) potential occurs at order e 2 of the perturbation theory. This 
distortion is more pronounced for larger values of E as shown in Figure 1. Note that this figure 
uses Eq. (59) which is free from the above-mentioned numerical error. 

I wish to thank Christiane Quesne for informing me of the above error. 
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